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The nonequilibrium real-time dynamics of electron coherence is explored in the quantum transport 
through the double-dot Aharonov-Bohm interferometers. We solve the exact master equation to find 
the exact quantum state of the device, from which the changes of the electron coherence through 
the magnetic flux in the nonequilibrium transport processes is obtained explicitly. We find that the 
relative phase between the two charge states of the double dot localizes to ? or — ? for all different 
magnetic flux. This nontrivial phase localization process can be manifested in the measurable 
occupation numbers. 
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Introduction. — Quantum coherence and quantum 
transport in mesoscopic electronic systems has attracted 
much attention due to recent achievements in quantum 
technology. Controlling quantum coherence in various 
quantum devices is essential to their functioning. Espe- 
cially in a ring structured Aharonov-Bohm (AB) inter- 
ferometer where electron interference can be tuned by an 
externally applied magnetic flux, a great amount of ef- 
forts has been made to investigate its coherence transport 
properties, mostly in the steady limit [H-flij. Alterna- 
tively we address here the real-time dynamics of building 
electron coherence and then the subsequent evolution of 
the intrinsic coherence in a two-terminal AB interferom- 
eter with two quantum dots (QDs). 

Coherence of electron transport through an AB inter- 
ferometer has long been characterized via conductance 
oscillation in magnetic fields [H-Q. Factors influencing 
AB oscillations include "which path detection" [H, 0, [|| , 
electron-electron interaction 0,0,0] and inelastic scatter- 
ing with phonons [loj . etc. Extracting the transmission 
phase from AB oscillations has always been the focus (see 
[lllll2| and references there in) since it may contain infor- 
mation about coherence of electron transport other than 
the interference by the AB phase. For a double dot AB 
interferometer, what most intuitively depicts the elec- 
tron coherence other than AB interference is the relative 
phase between the two dot charge states. This intrinsic 
phase is generally entangled with the AB phase and its 
fundamental dynamics has not been well understood so 
far. However, the dynamics of this intrinsic phase is im- 
portant for manipulating quantum coherence in the ap- 
plication of quantum information processing. Therefore 
in this Letter we shall directly study the nonequilibrium 
electron dynamics to examine the effects of the magnetic 
flux on the relative phase. 

To explore the intrinsic coherence dynamics, we solve 
exactly the time evolution of the electron charge states 
in the double dot AB interferometer. We found that 
although different values of magnetic flux will induce dif- 
ferent relative phases initially, they will eventually be 



localized to ir/2 or — ir/2. Remarkably, this phase lo- 
calization dynamics is reflected in the occupation num- 
ber rather than the transport current. For the latter, its 
magnetic flux dependence has been extensively employed 
in characterizing the coherence of electron transport. 

Reduced density matrix and exact master equation. — 
The system under consideration is described by the 
Hamiltonian consisting of three parts, H = rl s + %i ea d + 
rl c . rl s = ^ij a \ a j where i, j = 1,2 is the Hamilto- 
nian of the interferometer, with a i (a|) destroying (cre- 
ating) an electron on dot i. The lead Hamiltonian is 

^lcad = T,ak e a k c lk c c>k> where a = LoxR and c„ 
the electronic operators of the leads, and 
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jk 
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■H.c] (1) 



describes the coupling between the dots and the leads. In 
Eq. |T]) , <pj a are the AB phase factors due to the AB flux 

such that 4>il ~ (j)2L + 0i_r — 4>2R =<i>= 27T$/$ with 
$0 being the flux quantum. To single out the effects of 
the threading magnetic flux, we will not include explicitly 
the electron-electron interaction here. 

Coherence between the two electron charge states of 
the double dot is embedded in the reduced density matrix 
p of the double dot system, which can be obtained by 
tracing over the lead electrons from the total reduced 
density matrix p tot {t) = e" m ( t ~ t °Vtot(to)e m(t " to) . The 
(ex act) equation of motion for p has been derived in Ref. 
HI, U I- Assuming that initially, at t = to, the dots are 
uncorrelated with the leads then 

p = -i[H' s (t),p(t)] 

+ E {7« (t)(2a 3 p(t)al - a\a ]P {t) - p{t)a\ aj ) 

ij 

+ %(t){ajp(t)4 - a\p(t) aj - a\a jP {t) + p{t) aj a^ . 

(2) 

Here, H' s (t) = J^ij e ij(t) a l a j i s the renormalized Hamil- 
tonian of the double dot; e' and all other timc-dcpcndcnt 
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coefficients in Eq. ([2]) are given by eL-(f) — i%j(t) 



i[uu 1 ] ip ^Jt) 



v]^ , where the ma- 



trix functions u and v obey the following equations [1 



u{t) + iJ5u(r) + f dr'g(T - t')u{t') = , (3a) 

Jt 

v(r) +iEv(t) + ( dr'g(T -t')v(t') 

J t 



dr'g{T - r')u\t - r' + to) , (3b) 



to 



with the initial conditions u(to) = I,v(to) — 0. In 
Eq. (J3j) , is the 2x2 energy matrix of the double-dot, 
g and g are the temporal correlation functions resulted 
from the couplings to the leads via the spectral density, 

r m » = 2tt J2 V iah V jah e i ^-'-1'^6{u - e ak ) , (4) 

such that 0: g(r) = E,^, J F^aH an d 
9(r) = E Q= L,iJf^He^ T r Q M, where f a (cu) - 
(exp[(<J — /Lt Q )/fcBT] + 1) _1 is the initial Fermi distribu- 
tion for the electron reservoir (lead) a = L, R. 

Exact solution of the master equation. — To monitor 
the formation of coherence in the transport process, we 
prepare the double dot with the empty state 0} at t — to. 
We denote the singly occupied states by |1) and |2) re- 
ferring to occupation of the first dot and the second dot, 
respectively, as the two charge state basis. And the dou- 
ble occupancy is denoted by |3). From the master equa- 
tion ([2]) with the above preparation, the reduced density 
matrix becomes 



P(t) 



fp 00 (t) \ 

vu(t)-p3 3 (t) v 12 (t) 

v 2 i(t) v 22 (t) - p 33 (t) 

V o o o P33 (t)j 



(5) 



where the dynamics of the relative phase between the two 
charge states is given explicitly by the off-diagonal ma- 
trix element pnit) = Vi2(t), while poo(t) = dct[7 — v(t)] 
and ps3(t) = detv(t) determine the leakage effect and 
the double occupation, respectively, and pu(t) = vu(t) — 
P33(t) (i = 1,2) is the probability of each singly occupied 
charge state. It is obvious that the total probability is 
conserved: tr p(t) = 1. The full reduced density matrix is 
thus determined from v(t) which has a general solution 
from Eq. (l3b|). 



(0) 



V{t) = J — u(t,w)$^/ a (w)r B (w)ut(t,w) , 



where u(t,w) = J* dTe iu[t ~ T) u{t - r + t ) . It is then 
sufficient to solve the first line of Eq. ^ for a full under- 
standing of the intrinsic coherence dynamics. 



As being widely studied in the literature, with a de- 
generate double dot given by the Hamiltonian TL S = 
Ealcij and an equal coupling = = T/2 in the 
white band limit, we can well suppress the relaxation dy- 
namics and exclusively focus on the relative phase dy- 
namics. Also, as a convention, we choose the gauge 
4>il — 4>2L = 4*iR — 4>2R = </>/2. These considerations 



lead to Tl 



a ±i0/2- 

tain from Eq. ([3a]), by taking to = 0, 



r / 1 



One can easily ob- 



u(t) 



(cosh^ - a x S{(/))smhJ-^- 



(7) 



where T = T| cos(0/2)| and S(0) = Straight- 
forwardly one can also find from Eq. (|6]) 

v(t) = v (t)I + v x (t)a x + v y (t)a y + v z (t)a z , (8) 

where v a {t) = A + (t) + A_(t) and {v x (t), v y (t), v z (t)} = 
{S(<f>)(A + (t) - A_(t)), -ReB(t), S(<j))ImB{t)} with 

A t (t) = -(l±|cos||) / ^[f L {u) + f R {u>)]\u ± {t,u>)\\ 



B(t) = ^sm^ J ^[f L {w)-f R {w)]u + {t,w)u*_{t,u), 



(9) 



-E)t- 4(r±r , )t_ 



The notation cr, with i 



w ±(^ w ) = i (u -£)-i(r±r*) 
x 7 y,z denotes the Pauli matrices. 

Phase localization. — The electron coherence dynamics 
is embedded in the two charge state density matrix [the 
central block matrix in Eq. ([5])] which can be rewritten 
as 

P,(t) = l[l + r(t)-<r]- m{t)+ 2 m(t) I, (10) 

where r(t) = 2{v x (t),v y (t), v z (t)} is the polarization vec- 
tor of the two charge states, and poo(t) + P33(t) clearly 
shows the leakage effect. To exclusively study the co- 
herence between the two charge states, we shall apply a 
bias, pl — E = E — p R such that the two dots would 
be equally occupied, pn{t) = p22{t), i.e. r z (t) = be- 
cause Im(it + (i, uj)u*_{t 1 tu)) is antisymmetric in lu with re- 
spect to E. The polarization vector then is fully specified 
by r x (t) = 2Rep 12 {t) = 2v x (t) and r y (t) = ~2lmp 12 = 
2v y (t), which shows purely the dynamics of coherence be- 
tween the two charge states through the relative phase 
<p(t), defined explicitly by 

Pi2(i) = |Pia(*)|e ,v(t) = 5[r ie (t)-ir J ,(t)]. (11) 

Fig-fflplots the time evolutions pi 2 . At <j> = 2m-K where 
m is an arbitrary integer, Repi2 soon grows to stable 
values ±1/4 (with a time scale ~ and Impi2 keeps 
zero, which locks the relative phase <p to or tt. However, 
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when ^ 2m7r, Re/?i2 also grows to a maximal value 
(coherence building) with the same time scale and then 
decays to zero at a flux dependent rate T(l — | cos(0/2)|). 
Eventually all the different relative phases are localized 
to 7r/2 or —it/2. Fig. E] further visualizes this process. 
At short time t = 2/T, the direction of r(t) sweeps over 
the whole plane, which shows all kinds of relative phases 
between the two charge states induced by different values 
of the magnetic fluxes (— 2-k < 4> < 27r) [see Fig. EH a)], 
as a process of building coherence from the initial empty 
state. These different values (corresponding to different 
relative phases) then move toward the y axis as time goes 
on, except for the points <j> = (0, ±27r), sec Fig. E^b) and 
(c), and finally reside on it as the asymptotic limit shown 
in Fig. Eld) as the phase localization. 
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FIG. 1: The time evolution of pi2- (a) Repi? is plotted for <f> 
from to 2-7T since it is antisymmetric to (p = 7T. (b) Impi2 is 
plotted for (f> from — 7r to 7r and it is antisymmetric to cj> = 0. 
Here we take /xl — ^fl = 6F at temperature fcsT = T/5. 

The period of the two charge states in <f> is Air as a result 
of the intrinsic geometry of a two-level system. Though 
the relative phase is constant over the whole range of 
— 2-7T < (f) < 2tt, except for 4> = 0, the degree of interfer- 
ence 1/512 1 continuously changes with the flux as well as 
the bias. At zero bias, \pi2\ becomes zero in the steady 
limit, as a result of the statistical equilibrium. Therefore, 
the nonzero bias is crucial to manifest the above phase 
localization. We should also point out that this phase 
localization over the magnetic flux is gauge independent. 

To observe such process, we examine two measur- 
able quantities, the transport current passing through 
the double dot and the occupation number in the two 
dots. The former, usually derived from the nonequilib- 
rium Green function technique [ItJ , can also be obtained 



within our theory (see the explicit derivation in [14( and 
the result is consistent with [17]): 



Ht) = ~Re j ^(/lH - /«M){(1 + I cos ||)« + (t, u) 

+(1 — | cos ^|)u_(i, uj) — Tsin 2 ^u + (t,uj)ut(t, w)|. 

(12) 

The last term in the above current expression is indeed 
proportional to Im/9i2. The result is plotted in Fig.(3][a). 
As one can see, the current always oscillates smoothly in 
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FIG. 2: Distribution of (r x ,r y ) among different values of </> 
from — 2n to 2ir at different times. The parameters are the 
same as that used in Fig. [T] Note that changing the direction 
of the magnetic flux (f> — > — (f>, one flips the sign of r y while 
the sign of r x is flipped around <j> = ±7r as seen also in Fig. \T\ 



4> (showing the AB oscillation alone). The sharp transi- 
tions across </> = 2rmr and the invariance of the relative 
phase over the whole range of flux, as the main features 
of phase localization, do not be seen in the transport 
current. However, we find that the occupation num- 
ber, rii = pa + P33 with i = 1,2, which depends on 
Repi2 via hi(t) = -Tm(t) - V cos((j>/2)Rcp 12 (t) - %(t), 
does display the features of the phase localization dis- 
cussed above. The occupation number initially oscillates 
smoothly in <j), and in the steady limit, 



m(t -> oo) = - 



cos(</>/2)Re/5i2(£ — > oo). 



(13) 



As shown in Fig. |3jb), the phase localization to 7r/2 and 
— 7r/2 [as a result of Repi2(i — ► oo) = 0] are explicitly 
manifested by the invariance of rii(t — > oo) = 1/2 over 
the whole range of the flux, except for <fi = 2mir (m = 
0, ±1). At cf> = 2mn, sharp change of Repi2 to ±1/4 is 
also manifested as a sharp reduction of the occupation 
number from 1/2 to 1/4. 

Fig. 0(c) and (d) further reveals the underlying pic- 
ture of the phase localization. It shows that pn only 
contains the AB oscillation (oscillates in cf> smoothly), 
just like the transport current I does [see Fig. EJa) and 
(c)]. The asymptotic constant occupation 1/2 (mani- 
festing the phase localization) for any flux <\> ^ 2mn 
is contributed mainly from the double occupation P33. 
Furthermore, the sharp reduction of occupation number 
from 1/2 to 1/4 is indeed rooted in the vanishing double 
occupancy at <f> = 2rmr. As we know, strong inter-dot 
Coulomb repulsion may prohibit double occupancy re- 
gardless of the flux. Phase localization then might not 
be expected in this case. However, the parameters of the 
system can be well tuned to make the inter-dot Coulomb 
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repulsion become unimportant so that the phase local- 
ization can always be manifested. 
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FIG. 3: The dependencies of the current /, the occupation 
number n = ni = ni = pn+pas, single occupancy probability 
pn and double occupancy pss on <j) are examined at various 
times (the line styles are t = 0.6/r: pink ultrashort dashed, 
t = 1A/T: purple dash-dot, t = 2/Y: green dash-dot-dot, 
t = 6/F: green long dashed, t = 40/T: red solid and t — oo: 
blue short dashed). Both pn and P33 oscillate in flux in such 
a way that n is kept as a constant for (f> 7^ 2mir. 



The above discussion has centered on degenerate dou- 
ble dot. Practically there is no perfect degenerate double 
dot. However, we find that slightly splitting the degener- 
acy only slightly deviates the result of phase localization 
while the main property is still well preserved. As one 
can see from Fig. 0|a), r x = 2Repi2 is lifted a little bit 
from zero for the flux values near 2m7r. Similar signals on 
the occupation number can also be seen reflecting such 
a slight deviation for the phase localization, as shown in 
Fig. E|b). Comparing to the transport current, the elec- 
tron occupation in each dot is easier to be measured in 
experiments using, for example, SET or QPC [lij ]. this 
phenomenon can be well observed experimentally. 

Conclusion. — We have directly studied the dynamical 
effects of the magnetic flux on the relative phase between 
two electron charge states in a double dot AB interferom- 
eter by exactly solving the nonequilibrium electron dy- 
namics of the system. We found that the flux dependence 
of the relative phase characterizing the intrinsic elec- 
tron coherence lead to a phase localization through the 
nonequilibrium electron transport over the whole range 
of flux except for the points (f> = 2mir (m = 0, ±1 ). We 
also found that such distinguished dynamics of phase lo- 
calization is manifested in the occupation number which 
is expected to be measurable in experiments. 
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FIG. 4: The distribution of the polarization (r x ,r y ) and the 
occupation number n\ deviated away from the degenerate 
double dot, with energy level difference E\ — E2 = 2.5%r, 
are plotted in (a) and (b), respectively, at various times. The 
symbols and the corresponding time values used in (a) are the 
same as that used in Fig. [2] The different curves correspond- 
ing to various times in (b) are the same as that in Fig. G2b). 
712 is almost the same as ni since the deviation is small. 
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